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Abstract 

In this work, we develop a novel fault detection and isolation (FDI) scheme for discrete-time two- 
dimensional (2D) systems that are represented by the Fornasini-Marchesini model II (FMII). This is 
accomplished by generalizing the basic invariant subspaces including unobservable, conditioned invariant 
and unobservability subspaces of ID systems to 2D models. These extensions have been achieved and 
facilitated by representing a 2D model as an infinite dimensional (Inf-D) system on a Banach vector 
space, and by particularly constructing algorithms that compute these subspaces in a finite and known 
number of steps. By utilizing the introduced subspaces the FDI problem is formulated and necessary 
and sufficient conditions for its solvability are provided. Sufficient conditions for solvability of the 
FDI problem for 2D systems using both deadbeat and LMI filters are also developed. Moreover, the 
capabilities and advantages of our proposed approach are demonstrated by performing an analytical 
comparison with the currently available 2D geometric methods in the literature. Finally, numerical 
simulations corresponding to an approximation of a hyperbolic partial differential equation (PDF) system 
of a heat exchanger, that is mathematically represented as a 2D model, have also been provided. 
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I. Introduction 

Over the past few deeades, the problem of fault deteetion and isolation (FDI) of dynamieal 
systems has inereasingly reeeived larger interest and attention from the eontrol eommunity [[T| 
(c.f. to the referenees therein). The inereasing eomplexity of human-made maehines and deviees, 
sueh as gas turbine engines and ehemieal proeesses, has neeessitated the need for development of 
more eomplex and sophistieated FDI algorithms. Not with standing this requirement, development 
of FDI algorithms for systems that are governed by partial differential equations (PDF) has not 
been fully addressed and reeeived extensive attention in the literature. One approaeh to investigate 
the FDI problem of PDFs relies on obtaining an approximate model of the system. First, the 
PDE system is approximated by a simple model (sueh as an ordinary differential equation (ODE) 
model), and then suffieient eonditions for solvability of the EDI problem are derived based on 
this approximate model. 

It is well-known that parabolic PDE systems ean be approximated by ODE representations. 
These systems ean be approximated through applieation of finite element methods where suffi¬ 
eient eonditions ean then be derived by using singular perturbation theory Q . The FDI problem 
of parabolie PDFs has been addressed by using the eorresponding approximate models in 0, 0. 
On the other hand, by diseretizing through spatial eoordinates, one ean approximate hyperbolic 
PDE systems by ODE models. However, unlike parabolie PDE systems, one eannot apply model 
deeomposition, order reduction and singular perturbation theory to hyperbolic PDE systems 
0. Moreover, the order of the resulting approximate ODE systems can be dramatically high. 
Therefore, researchers have investigated hyperbolic PDE systems by using other formal methods 
such as the theory of semigroups 0 and backstepping methods 0. 

As shown in 0, a single hyperbolic PDE system can be approximated by using two-dimensional 
(2D) systems. In 0, we have shown that this approximation is also applicable to a system of 


hyperbolic PDFs. Moreover, as shown in [10|, parabolic PDE systems can be approximated by 


three-dimensional (3D) Fornasini and Marchesini model II (FMII) representations [ 111-[ 14|. We 
have shown in the Appendix that it is straightforward to extend the results of our work to 3D 
systems as well (these results are not included here due to space limitations). Consequently, our 
proposed methodology in this work can also be applied to parabolic PDE systems. There are 
only a few results on FDI of 2D systems in the literature, such as those by using dead-beat 
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observers p3| and parity equations p^ , p7| . On the other hand, the FDI of hyperbolic PDE 
systems (by using even an approximate model) has not yet been adequately addressed in the 
literature. 

It should be noted that 2D system theory has other applications in the control field. For 
example, a class of discrete-time linear repetitive processes can be modeled by 2D systems. 
These processes play important roles in tracking control and robotics, where the controlled 


system is required to perform a periodic task with high precision (refer to [ 18| for more details 
on repetitive systems). One of the main approaches to control linear repetitive processes is the 
iterative learning control (ILC) p^. Since the ILC problem can be formulated as a control 


problem in the 2D system theory [19|, 2D systems have been increasingly applied to spatio- 
temporal and repetitive process control problems in the literature. 

The FDI of ID linear, time-invariant systems has been extensively investigated during the past 
few decades [|T| (and the references therein). The geometric approach p0|-p2| has provided 
a valuable tool for investigating the FDI problem of a large class of dynamical systems such 
as linear time-invariant 


|, Markovian jump systems p3[ , time-delay systems p4| |, linear 
impulsive systems [ |^ , and parabolic PDE systems 0. Moreover, the geometric approach is 
also extended to affine nonlinear systems in p6| . Furthermore, hybrid geometric FDI approaches 


for linear and nonlinear systems have also been provided in [271 and [^, where a set of residual 
generators are equipped with a discrete event-based system fault diagnoser to solve the FDI 
problem. 

Motivated by the above discussion, in this paper we investigate the FDI problem of 2D systems 
and apply the results to a 2D approximate model of a hyperbolic PDE system. As stated earlier, a 
hyperbolic PDE system can be approximated by ID systems, where the order of the approximate 
ID system will significantly increase by decreasing the gridding size Q. We also provide an 
example where faults in the ID approximate model are not isolable, whereas one can detect and 
isolate faults by applying the 2D approximate model and by using our proposed methodology. 

Recently, the geometric theory of 2D systems has attracted much interest, where basic concepts 
such as conditioned invariant and controllable subspaces are studied in detail for the Fomasini 


and Marchesini model I (FMI) [29|, |30|. The hybrid 2D systems have also been investigated 


from the geometric point of view in [311. The geometric FDI approach for 2D systems, for the 
first time, was addressed in Q, where invariant subspaces of the Roesser model are defined and 
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the FDI problem is formulated based on these subspaees. 

Two-dimensional (2D) systems have been extensively investigated from a system theory point 
of view I[TT|-||T4|. Partieularly, system theory eoneepts sueh as stability p3| |, [ [^ , eontrollability 
p3[ , observability pT| , and state reeonstruetion p4| have been investigated in the literature. 
However, due to eomplexity of 2D systems, unlike one-dimensional (ID) systems, there are 
various definitions that are introdueed for eontrollability and observability properties. Not sur¬ 
prisingly, the duality between the observability and eontrollability does not hold in 2D systems. 
In this paper, we investigate observability of 2D systems from a new geometrie point of view 
whieh has its roots in system theory of infinite dimensional (Inf-D) systems. 

Compared to the results reported in Q and p5| , we have speeifie generalization and novel 
eontributions in this work. We first investigate the Fomasini and Marehesini model II (FMII) as 
an Inf-D system that allows us to deal with Inf-D subspaees (albeit with a finite dimensional (Fin- 
D) representation). The invariant subspaees and the eorresponding algorithms that are introdueed 
m iivii for the Roesser model are then generalized to the 2D FMII systems. It is worth noting 


that although the introdueed subspaees in this work are Inf-D, the eorresponding algorithms for 
eonstrueting the subspaees eonverge in a finite and known number of steps. 

In addition, in [[^ only suffieient eonditions for solvability of the FDI problem of the Roesser 


model were provided. As shown in |35|, the invarianee property of an unobservable subspaee is 
a generie property of 2D systems. We first derive a single neeessary and suffieient eondition for 
deteetability and isolability of faults (that is formally introdueed in the next seetion). It is shown 
that this eondition also neeessary for solvability of the FDI. Furthermore, in 0 by utilizing 
the existenee of an LMI-based observer only suffieient eonditions for solvability of the Roesser 
model FDI problem were provided. In other words, the proeedure to design the observer gains 
is not provided in Q, [ [35| . However, in our paper, we derive both neeessary and suffieient 
eonditions, where suffieient eonditions are based on (a) an ordinary, (b) a delayed deadbeat, and 
(e) an LMI-based 2D Luenberger filters. Moreover, we develop a proeedure to design LMI-based 
filter gains. 

It must be noted that reeently related work has appeared in [[3^ and [|37|. These two papers 


investigate the FDI problem of three-dimensional (3D) FMII models. Although, a geometrie FDI 
methodology is also developed in [36|, our work is distinet and unique from [ [3^ in three main 
perspeetives: 
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1) The approach proposed in is based on the results of pT] |, whereas our approach is 

(the results in 


based on the generalization of the results of pT| as reported in 
are more general than [211) for 2D models (for more information refer to Remark [^. 


2) In [361, necessary and sufficient conditions for solvability of the FDI problem were derived 
for a subclass of detection filters where it was assumed that the output map of the detection 
filters and that of the system are identical. However, in our work here we consider a 
general class of detection filters for the residual generation and relax this condition. 


3) As shown in Section IV-B the observability property of the 2D model is a fundamental 


requirement and assumption in [361 (although it is stated in [36| that this assumption was 
made for simplicity of their presentation). However, our proposed solution does not require 
this condition and assumption, and consequently our approach leads to a less restrictive 
solution. 

Another approach that was developed in the literature [ |T5| , [ [T7| has its roots in 2D deadbeat 
observers 


|. In [ [T5| [, the FDI problem is investigated by using polynomial matrices and 
unknown input deadbeat observers, where the right zero primeness of the 2D Popov-Belevitch- 
Hautus (PBH) matrix (this is reviewed in Subsection [II-D[ ) is a necessary condition for solvability 
of the FDI problem. In [ [T7| , this condition was relaxed and necessary and sufficient conditions 
that are based on an extended parity equation approach were obtained. 


To provide a fair and comprehensive comparison with the currently available result in [36[, in 
this work it is shown that on one hand the solvability of the FDI problem by using the method 


in [36[ is also sufficient to accomplish the FDI task by using our proposed approach. On the 


other hand, there are certain 2D systems that are not solvable by the approach in [36[, however 
our proposed approach can both detect and isolate the faults. Also, by comparing our proposed 
results with the algebraic-based methods in [ [T5| , [ [T7| [, where they can solve and derive necessary 
and sufficient conditions for solvability of the FDI problem, we will highlight and emphasize 
two important considerations as follows: 

1) The algebraic methods, in contrast to our geometric approach, need a closed-form and 
analytical solution to certain polynomial matrix equations (in two variables). However, 
our proposed approach is derived and solved by using commonly available and relatively 
straightforward numerical methods. 
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2) As shown subsequently in Seetion IV-B there are eertain examples, where the nee- 
essary eonditions in [15|, 10 are not satisfied. However, our proposed approaeh ean 
both deteet and isolate these fault seenarios. 

To summarize, the main eontributions of this paper ean be highlighted as follows: 

1) By reformulating 2D models as Inf-D systems, the invarianee property of the unobservable 
subspaee is investigated (this is provided in Seetion |I^, where an Inf-D unobservable 
subspaee is also introdueed. This result enables one to formally address the solvability of 
the FDI problem without restrietion on the initial eonditions (unlike in Q, [ [^ where 
a restrietive assumption that the unobservable subspaee is Fin-D is imposed). Sinee 2D 
systems are Inf-D dynamieal system, our proposed Inf-D representation and framework 
enables one to address the FDI problem in its most general seenario than that in Q and 


m. 

2) Two important Inf-D invariant subspaees (namely, the eonditioned invariant and the unob¬ 
servability) are introdueed for the FMII models. Although, these subspaees are Inf-D, we 
provide explieit algorithms that ean be invoked to eompute these subspaees in a finite and 
known number of steps. 

3) The FDI problem of 2D systems is formulated in terms of the above introdueed invariant 
subspaees, and neeessary and suffieient eonditions for its solvability are derived and 
formally analyzed. 

4) A novel proeedure is developed for designing an observer (also known as a deteetion filter) 
by utilizing the linear matrix inequalities (LMI) teehnique. 

5) Three sets of suffieient eonditions for solvability of the FDI problem by utilizing the 
ordinary, the delayed deadbeat, as well as our proposed LMI-based observers (deteetion 
filters) are also provided. 


6) Analytieal eomparisons between our proposed approaeh and the one in [36| are presented. 


We show that the suffieient eonditions in [361 are also suffieient for solvability of the 
FDI problem by using our approaeh. However, an example is provided that shows our 
method ean both deteet and isolate faults, whereas the approaeh in eannot be used. In 


other words, it is shown that if the method in the above literature ean aeeomplish the FDI 
task, our proposed approaeh ean also aeeomplish this task. However, if our seheme eannot 
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achieve the fault deteetion and isolation goals for a given system, then it is guaranteed 
that the sehemes in [|^ eannot also aehieve these goals. Moreover, there are 2D systems 


where our approaeh ean aehieve the FDI objeetives whereas the results in the literature 
eannot solve the FDI problem. 

7) Our proposed methodology and strategy is applied to an important applieation area of a 
heat exehanger (a hyperbolie PDF system), where it is shown that one ean simultaneously 
deteet and isolate two different faults namely, the leakage and the fouling faults. 

The remainder of the paper is organized as follows. The preliminary results ineluding the Inf-D 
representation, the FDI problem formulation, the 2D deadbeat observer and the 2D Luenberger 
observers (deteetion filters) are presented in Seetion The unobservable subspaees of the 


FMII 2D model are introdueed in Seetion III The geometrie property of these subspaees 
and the invariant eoneept of the FMII model are also presented in Seetion III In Seetion 
neeessary and suffieient eonditions for solvability of the FDI problem are derived and developed. 
Analytieal eomparisons between our proposed approaeh and the available geometric methods 


in the literature, namely [361 and [37| are provided in this seetion. Furthermore, numerieal 
eomparisons with both geometric and algebraic methods in [15|, 10, [ |^ , [ [37| are presented 
in this seetion. Simulation results for the FDI problem of a heat exehanger that is expressed as 
a PDF system are eondueted in Seetion |Vj Finally, Seetion VI eoneludes the paper and provides 
suggestions for future work. 

Notation: In this work, are used to denote subspaees. For a given veetor L, the 

subspaee span{L} is denoted by The inverse image of the subspaee 7^ with respeet 

A 0 

is denoted by 


to the operator A is denoted by A . The bloek diagonal matrix 


0 B 


diag(74, B). The real, eomplex, integer and positive integer numbers are denoted by M, C, Z and 
N, respeetively. N denotes the set N U {0}. In this paper, we deal with infinite dimensional (Inf- 
D) subspaees and veetors. An Inf-D veetor is designated by the bold letters x, y, • ■ ■. The Inf-D 
subspaee ■ • •©7^©T'©- ■ ■ is denoted by ©7^, where Y C M”. Let x = (■ ■ • , , xj, ■ ■ C 

©7^ and |x|c 
©IL and |x 

spaee. Let i, j e Z U {—cxd, cxd} and J > The Inf-D veetor x^ E ©7^ is expressed as x^ = 
[• ■ ■ , 0, 0, xj, • • • , xj, 0, 0, • • ■ where xe E Y for alH < £ < j, and assoeiated with we 


= sup|xi|, where x* E 7^. The veetor spaee 7^ = is defined as {x|x G 

ieZ 

< cxd}. It ean be shown that 7^ is a Banaeh (but not neeessarily Hilbert) 
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simply use x. The other notations are provided within the text of the paper as appropriate. 

II. Preliminary Results 

In this seetion, we first review 2D systems and their various representational models. Sub¬ 
sequently, a 2D system is expressed as an infinite dimensional (Inf-D) system that allows one 
to geometrically analyze the unobservable subspaces (this is to be defined and specified in the 
next section). The FDI problem is also formulated in this section. Moreover, we review the 
2D Popov-Belevitch-Hautus (PBH) matrix and 2D deadbeat observers in this section. Finally, 
an LMI-based approach is introduced to design a 2D Luenberger observer (also known as a 
detection filter) for 2D systems. 


A. Discrete-Time 2D Systems 


2D models can be used for representing a large class of problems such as approximating hyper¬ 
bolic PDF systems @, @, image processing and digital filtering |j^. System theory concepts 
such as observability, controllability and feedback stabilization have also been investigated in the 
literature for 2D systems Q, pT|, [14|, [30|, [311. There are various models that are adopted in 


the literature for 2D systems including the Rosser model |39|, the Fornasini-Marichesini model 
I (FMI) and model II (FMII) [ [TT| , [ [T4l . The FMI can be formulated as a Roesser model and 
the Roesser model is a special case of the FMII model [ [TT| . In this work, we consider and 
concentrate on the FMII model, and consequently our results are also derived for this general 
class of 2D systems. 

Consider the following FMII model [ [T4| , 

x{i + IJ + 1) = Aix{iJ -f 1) + A 2 x{i -f 1, j) -f Biu{iJ + 1) -f B 2 u{i + 1, j) 


^ Lkfkihj + 1) + ^ Llfk{i + 1, j), 


( 1 ) 


k=l 


k=l 


yihj) = Cx{i,j), i,j e z, 

where x G M”, u G M™, and i/ G denote the state, input and output vectors, respectively. The 
fault signals and the corresponding fault signatures are designated hyfk, L], and L^, respectively. 
Also, p denotes the number of faults in the system. Since in this work all the introduced invariant 
subspaces are based on the operators Ai, A 2 and C, we designate the system ([T]) by the triple 
(C',Ai,A2). 
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Remark 1. Note that system ([T]) represents and captures the presence of both actuator and 
component faults. To represent sensor faults, one can augment the sensor dynamics and model 
the sensor faults as actuator faults in the augmented system (for a complete discussion on this 
issue refer to l \20^ - Chapters 3 and 4). Also, it should be pointed out that the fault signal fk 
affects the system through two different fault signatures L\ and L\. An alternative fault model 
could have been expressed according to the following representation, 

p 

x{i + l,j + 1) =Aix{i,j + 1) + A 2 x{i + 1, j) + + 1) + B 2 u{i + l,j) + '^Lkgk{i,j), 

k=l 

y{i,j) =Cx{i,j). (2) 

Model ([T]) is more general than the one given by equation This is due to the fact that by 
denoting fk(i + 1, j) = gk{i,j) for all k = 1, ■ ■ ■ ,p, one can represent the model © as in the 
model O- ■ 


Let us now consider the Roesser model [|39| which is expressed as 


r(i + l,j) 


2I11 

Ai2 


i 

+ 

Bii 






s{i,i + 1 ) 


A 21 

A 22 




B 21 


u 


( lj ) N^LkfSd), 


k=l 


y{hj) =c 


s{hj) 


(3) 


and where 


1 T 


T T 

r s 


represents the state, and the variables u, y, fk and Lk are defined 


as in equation ([T]). By defining 


X = 


r 

5 Ai — 

1 

(N 

, ^ 

5 A 2 — 

1 

0 

0 

,5l = 

Bii 

,i ?2 = 

0 

s 


1 

0 

0 


A 2 I A 22 


0 


B 21 


(4) 


one can formulate the Roesser model @ as in equation ([T]). In this paper, we assume that Ai 
and A 2 in model ([T]) are not necessarily commutative (i.e. ^ 1^42 7 ^ A 2 A 1 ), and hence, the results 
that are subsequently developed can also be applied to the Roesser model It should be noted 


that the commutativity of Ai and A 2 is a strong condition that renders the results in [311 (where 
it is assumed that Ai and A 2 are commutative) not applicable to the system 

In this work, we will investigate and develop FDI strategies for the model ([T]). It is assumed 
that Ai and A 2 in model ([T]) are not necessarily commutative (i.e. A 1 A 2 7 ^ A 2 A 1 ), and hence. 
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the results that are subsequently developed ean also be applied to the Roesser model. It should 
be emphasized that the eommutativity of Ai and A 2 is a strong eondition that renders the results 


in [311 (where Ai and A 2 are assumed to eommutate) not applieable to Roesser systems. 


B. Infinite Dimensional (Inf-D ) Representation 

In this subseetion, we reformulate the 2D model ([T]) as an Inf-D system that will be used to 


derive the invarianee property of unobservable subspaees (for details refer to Seetion III-A l. 

Consider the fault free system ([T]), that is with /^ = 0, and with zero input (we are mainly 
interested in the unobservable subspaees and do not need to be eoneerned with the eontrol inputs 
in the FDI problem). By eonsidering x(/c) = (■••, x{—l + k, I)"*", x{k, 0)"’", x{l + k, —1)"’", •••)"'"€ 
^ it ean be shown that under the above eonditions the system ean be represented as, 

x(A; -f 1) = Ayi{k)^ k e N 

(5) 

y{k) = Cx(fc), 

where x(A;) g y(^) = (•••,2/(-l + k, 1)'^, y{k, 0)'^, y{l + k, e 

and A is an Inf-D matrix with Ai and A 2 as diagonal and upper diagonal bloeks, respeetively, 
with the remaining elements set to zero, and C = diag(- ■ ■ , C, C, • ■ ■). In other words, we have. 




■ ■ ■ 0 Ai A2 0 ■ ■ ■ 

,c = 

0 

0 

■ ■ ■ 0 0 Ai A2 ■■■ 

... 0 0 c ■■■ 





( 6 ) 


Note that sinee we invoke an Inf-D representation to investigate an unobservable subspaee, 
and where this subspaee is defined by only A and C, therefore for sake of presentation simplieity, 
an Inf-D system is used that has no fault and zero input. 

There are various formulations for the initial eonditions of the FMII model Q. These are 


based on the separation set that is introdueed in [40|. There are two separation sets that are 
eommonly used in the literature. In the first formulation the initial eonditions are denoted by 
x(0) = (••• , x(—1,1)"’", a;(0, 0)"’",x(l, —1)"’", ■ ■ ■)"’" G (this is eompatible with the 


model @). The seeond formulation is expressed as a;(i,0) = hfii) and x(0,j) = h2{j), where 
hi{i), h2{j) G MA and i, j G N 0 . The seeond formulation is more eompatible with applieations 
(partieularly, in ease that the system ([T]) is an approximate model of a PDF system - refer to 
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Section |V]). It will be shown subsequently that since we derive the conditions based on invariant 
unobservable subspace (this is formally defined in the next section), our proposed methodology is 
applicable to both initial condition formulations. In other words, we use the Inf-D representation 
to only show the results and evaluate the developed algorithms. However, to apply our results 
there is no need to deal with Inf-D systems and subspaces, and therefore, one can apply our 
proposed methods to 2D systems corresponding to both initial condition formulations. 

We start with the first formulation of the initial conditions subject to the boundedness assump¬ 


tion (this is, x(0) G ^M”). However, as shown in Section III-B our proposed results also hold 
for the second initial condition formulation. 

As stated in the Notation section, it can be shown that X defined for equation @ is an Inf-D 
Banach space. The system theory corresponding to Inf-D systems is more significantly challeng¬ 
ing than Fin-D system theory (ID systems) (refer to [[4T]|). However, as shown subsequently, the 
operator A is bounded and consequently, one can readily extend the result of ID systems to the 
system ([^ [[^, [|4T|. Let us first define the notion of bounded operators. 


Definition 1. Consider the operator A \ Xi ^ X 2 , where Xi and X 2 are Banach vector 
spaces with the norms | ■ |i and 1-12, respectively. The operator A is bounded if there exists a 
real number G such that |.Ax |2 < G|x|i for all x G Xi. 


Lemma 1. The operator A as defined in the Inf-D system @ is bounded. 


Proof: Let G = 2max(|Ai|, IA 2 I), where \ Ai\ denotes the norm of Aj and x = {xk)k& ^ A". 

It follows readily that |Ax|oo = suplAiXfc -f ^ 2 X^+ 1 ! < supGmaxdxfcl, Ix^+il) = Gsuplx/d. 

fcez fcez fcez 

Therefore, |Ax|oo < G|x|oo. This completes the proof of the lemma. ■ 


The above lemma enables one to now formulate the unobservable subspace of the 2D system 


([T]) in a geometric framework (for details refer to Section III) based on the operator A (and 


consequently, in terms of Ai and A 2 ). 


Remark 2. Although, in / [4^ all the results such as the controlled invariant subspaces are 
presented on the developed approach in / |?6] / has its roots in the theory of systems over 
rings. In this paper, we propose an alternative approach that is based on Inf-D systems that 
are defined on a Banach vector space. Similar to /[?^, our proposed methodology including the 
algorithms and the conditions for solvability of the FDI problem can also be addressed in a 


September 20, 2016 


DRAFT 






12 


Fin-D scheme. However, as shown in the literature the duality property does not hold for 2D 
systems /[74|/, pS]/. Therefore, by simply invoking duality the results of this paper cannot be 


derived from those in 146]. 


C. The FDI Problem of 2D FMll Model 

In this subsection, we formulate the FDI problem for the 2D system ([T]). In this paper, without 
loss of any generality, it is assumed that the system ([T]) is subject to two faults, and therefore 
we construct two residuals such that each one is sensitive to only one fault and is decoupled 
from the other. 

More precisely, consider the faulty FMII model ([T]). The solution to the FDI problem of the 
2D FMII system can be stated as that of generating two residuals rk{i,j),k G {1, 2} such that, 

Vu, /2 and /i = 0 then ri —)■ 0, 

and if /i 7 ^ 0 then ri 7 ^ 0, (7a) 

Vm, fi and /2 = 0 then r 2 —)■ 0, 

and if /2 7 ^ 0 then r 2 7 ^ 0. (7b) 

The above residuals are to be constructed by employing fault detection filters. For the 2D 
system 0 , we consider the following FMII-based fault detection filter, 

<4{i + 1,7 + 1) = Fiui{i,j + 1) + F 2 U){i + l,j) + + 1) + K 2 u{i + 1,7) + Eiy{i,j + 1) + E 2 y{i + 1,7), 

rk{i,j) = Muj{i,j) - Hy(i,j), (8) 

where uj{i,j) denotes the state of the filter and is used to define the residual signal rk{i,j). The 

solution to the FDI problem is now reduced to that of selecting the filter gains Fi, F 2 , Ki, K 2 , 

El, E 2 , M and H corresponding to the filter 0. 


Remark 3. The detection filter 0 can be selected as a full-order (H = I) or as a partial- 


order (kerH f 0) 2D Luenberger observer. As shown subsequently in Section Ty_ this level of 
generality allows one to analytically compare our proposed methods with those results reported 
in ( |36] /. ■ 


Remark 4. In this paper, we investigate the FDI problem by employing two main steps, namely 
(i) decoupling the faults, and (ii) designing filter gains for each fault. The first step for decoupling 
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/i addresses the existence of three maps Di, D 2 and H, such that the fault /2 signatures L\ and 
Li are members of the unobservable subspace (defined in the next section) of the system (HkC, 
Ai + DiC, A 2 + D 2 C). The same terminology is used to decouple / 2 . Moreover, the second step 
is mainly concerned with existence of the filter ([^ such that stability of the error dynamics is 
guaranteed. In this paper, if the first step is solvable for the fault fi we say that fi is detectable 
and isolable. Finally, it is stated that there is a solution to the FDI problem if for all the fault 
signals fi both steps above are solvable. 


D. Deadbeat Observers 


In Section IV, necessary and sufficient conditions for solvability of the FDI problem are 
derived. We provide sufficient conditions for accomplishing the FDI task by using a delayed 
deadbeat detection filter and an ordinary (i.e., without a delay) deadbeat observer (refer to the 
subsequent Corollaries and |^. Towards these end, in this section we formally define a (delayed) 


deadbeat filter. For a comprehensive discussion on 2D deadbeat observers refer to [17|, [381. 

Consider the system ([1]) under the fault free situation. A (delayed) deadbeat observer is 
constructed according to, 

zii + + 1) =Fiz{i,j + 1) + F2z{i + l,j) + Kiu{i,j + 1) + K2u{i + 1, j) + Eiy{i,j + 1) + E2y{i + l,j), 

=Liz{i,j) + L2yii,j), (9) 


where u and y are the input and output signals as defined in the system ([T]). Note that since 
the output is assumed to not be directly affected by the input signal, x{i,j) is only a linear 
combination of z{i,j) and y{i,j). If there exists a number Nq such that x{i,j) = x{i,j) for 
all i + j > Nq, the filter Q is designated as an ordinary (without a delay) deadbeat observer. 
On the other hand, if there exist non-negative integers rii and n 2 such that ni + n 2 > 0 and 
x{i,j) = x{i — ni, j — 712), the observer @ is designated as a delayed deadbeat observer. The 
necessary and sufficient conditions for existence of a (delayed) deadbeat observer are specified 
in the following theorem in).® 


Theorem 1. Consider the 2D system 0 under a fault free situation and the following 2D 
Popov-Belevitch-Hautus (PBH) matrix. 


PBH{zi,Z2) 


I — ziAi — Z2A2 
C 


( 10 ) 
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where Ai, A 2 and C are defined as in equation Q. and Zi,Z 2 G C. Then, 

(i) there exists a delayed deadbeat observer if and only if iank{PBH{zi, Z 2 )) = n for all 
Z\,Z 2 G C — {0} (that is, PBH{zi, Z 2 ) is right monomic) / |77| / (Section 3) and 

(ii) there exists an ordinary deadbeat observer if and only if iai]k(PBH(zi,Z 2 )) = n for all 


Zi,Z 2 G C (that is, PBH{zi, Z 2 ) is right zero prime) [38], 


We use the above theorem in Section |IV] to show the existence of deadbeat filters to derive 
sufficient conditions for solvability of the FDI problem (the subsequent Corollaries and [^. 


E. LMI-based Observer (Detection Filter) Design 

As shown in p^, design of a deadbeat observer Q requires that one works with polynomial 


matrices (this is not always a straightforward process). In this subsection, we address the design 
process for the FMII system observer, or the detection filter gains, by using linear matrix 
inequalities (LMI). These results will be used to explicitly design a 2D Luenberger detection 
filter (that can also be formulated as in equation @) subsequently in Section]^ for the purpose 
of accomplishing the solution to the FDI problem. 

In order to show the asymptotic stability of the state estimation error dynamics, one needs to 
apply the following stability lemmas. 


Lemma 2. The 2D FMII system ([T]) (under the fault free situation) is asymptotically stable 
if there exist two symmetric positive definite matrices R\,R 2 G such that. 


A, = A^(Ri + R2)A - i? < 0, 

where A = [Ai Af\ and R = diag(i?i, R 2 ). 


( 11 ) 


Lemma 3. [42] Consider the FMI condition <I) + P^IsAQ + KP < 0, where <I> G 
P G and Q G There exists a matrix A G satisfying the previous FMI condition 

if and only if Wp^Wp < 0 and Wj^Wq < 0, where the columns of Wp and Wq are bases of 
the ker P and ker Q, respectively. ■ 

Now consider the 2D system ([T]) under the fault free situation and the corresponding state 
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estimation observer as given by, 

+ Iji + 1) =(^i + DoiC)x{i, j + 1) + {A 2 + Do2C)x{i + l,j) + Biu{i,j + 1) 

- Doiy{i,j + 1) - Do2y{i + l,j) + B2u{i + 1, j), 

Vihj) =Cx{i,j). (12) 

It follows readily that the state estimation error dynamies, as defined by e{i,j) = x{i,j) —x{i,j), 
is governed by, 


+ I 5 J + 1) —(^1 + DoiC)e{i, j + 1) + {A 2 + Do2C)e{i + 1, j). (13) 

The following theorem and eorollary provide an LMI-based eondition for existenee of the state 
estimation observer gains Doi and D 02 sueh that the error dynamies ( [T3| ) is asymptotieally stable. 

Theorem 2. Consider the 2D system ([T]) under the fault free situation. There exist two maps 
Doi , Do 2 : —)■ and two symmetric positive definite matrices Ri and R 2 such that the 

LMI ( [TT] ) is satisfied for Ai + DoiC and A 2 + D 02 C if and only if Ri and R 2 satisfy the LMI 
condition Wj^AcWc^ < 0, where = dieig{Wc,Wc) and the columns ofWc € are 

the basis o/ker(7. 


Proof: Note that without loss of any generality, it is assumed that C is full row rank, and 

-(i?l -r rL2) 


n > q that is equivalent to partial state measurement. Let $ = 


R2)-^ a 




-R 


the Sehur eomplement lemma, we have WjAcWc^ < 0 if and only if. 


By using 


— (/?! + R 2 ) ^ 



I 

0 

$ 

I 

0 


-WlRW,,_ 


0 



0 



< 0 . 


(14) 


It follows that 


2nxn ^2nx2n 


$ 


Onx 2n 
X 2n 


= —R < 0 if and only if i? > 0 (or > 0 and 


R 2 > 0). By defining P = 
and using Lemma the LM 


2qxn 


c. 


and Q = 


^nxn ^nx2n 


, where Cd = diag(C, C), 


A = 


Doi Do2 


$ 


eondition ([T4|) is satisfied if and only if there exits a matrix 


G sueh that. 


0nx2q 

A T 



^nxn 




{Ri + R 2 ) ^ G 



^nxn 0nx2n 

+ 


A 

02(jxn Cd 

= 


1- 

1 _ 




^2nxn 




1 

q 


< 0 , 


(15) 
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where G = 


+ DoiC A2 + D02C 


. Again, by using the Sehur eomplement lemma, we have 


G^{Ri + R2)G - R < 0. This eompletes the proof of the theorem. 

An important eorollary to the above theorem and Lemma |^ean be stated as follows. 


Corollary 1. Consider the 2D system o under the fault free situation and the state estimation 
observer ( [T^ . If there are two symmetric positive definite matrices Ri and R 2 satisfying the LMI 
condition W'^^AcWc^ < 0, then there exists two maps Doi and D 02 such that the error dynamics 
0 is asymptotically stable. 


Proof: Follow direetly from Theorem and Lemma and the details are omitted for sake 
of brevity. ■ 


Remark 5. Note that by solving the LMI condition Wj^AfWc^ < 0, one can obtain symmetric 
positive definite matrices Ri and R 2 . Hence, the state estimation observer gains Doi and D 02 
are computed by solving the equation (which is an LMI condition in terms of the gains Doi 
and Do 2 ). Therefore, Corollary not only provides sufficient conditions for existence of a state 
estimation observer, but also provides an approach for computing the observer gains Doi cind 
Do2- ■ 


The results of this seetion will now be used subsequently in Seetions III-A and IV to address 
the unobservable subspaee of the system ([T]) as well as to provide suffieient eonditions for 
solvability of the FDI problem, respeetively. 


III. Invariant Subspaces for 2D FMII Models 

As deseribed earlier, 2D systems ean be represented as Inf-D systems (i.e. the initial eondition 
is a veetor of an Inf-D subspaee). In this seetion, we first use the Inf-D representation Q to 
formally define and eonstruet an unobservable subspaee. Next, we define a subspaee of the 
unobservable subspaee (this we ealled as an invariant unobservable subspaee) of the 2D system 
0 that ean be represented as an infinite sum of the same finite dimensional subspaees. Therefore, 
one ean eompute the invariant unobservable subspaee (that is, the Inf-D subspaee) in a finite 
number of steps. Also, it is shown that the invariant unobservable subspaee enjoys an important 
geometrie property that is erueial for solving the FDI problem. 
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A. Unobservable Subspace 

As described in the previous section, 2D systems can be represented as Inf-D systems. In 
this subsection, the Inf-D representation Q is utilized to formally define and construct an 
unobservable subspace. 

The unobservable subspace of the system Q (and consequently of the system ([T])) is defined 
as, 

OO 

kerCAl^ ( 16 ) 

fc =0 

where A and C are defined as in equation @. Note that we define the above unobservable 
subspace by following along the steps in [ [4^ , the results in [ |4T| (Chapter I), and the fact that 
the operator A in equation (|^ is bounded (refer to Lemma [^. 

One of the main difficulties in geometric analysis of Inf-D systems is the convergence of any 
developed algorithm that involves computation of certain set of subspaces in a finite number of 
steps. For example, consider the unobservable subspace ( [T^ . In Fin-D systems, the algorithm for 
computing the unobservable subspace converges in a finite number of steps [ [44| . Moreover, one 
is generally interested in investigating the FMII models in a Fin-D representation ([T]). Motivated 
by the above, below two important subspaces of A^g that are denoted by and A^s,oo are 
introduced. The subspaces and A's^oo can be computed in a finite number of steps and also 
allows one to derive necessary and sufficient conditions for solvability of the FDI problem. 

{ Mo ; i = j = 0 

, 

0 ; otherwise 

where Uq G MA. One can show that the state solution of the model ([IJ) under the fault free 
situation is given by pA] |, 




(17) 


(18) 


where the matrices A67 )’s and A^’-^^’s are defined by the following recursive expressions, 

=0 if i or j < 0, 

Agj) _ ^6-1,A(°’°) = /. 

Based on the solution that is given by equation ( [TT] ), and considering that mq = 0, a finite 
observability matrix (given that its null space is a finite dimensional subspace) can be defined 
as follows. 


O = 


CT,(CAi)T,(C'A2)T,-- - ,(C'A(*’^))T, 


(19) 
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Let ^ = ker O = Hi j>o (kerCA^*’-^)). Since dim(^) < n < oo, we designate ^ as ihe finite 
unobservable subspace of the system ([T]). Also, recall from the 2D Cayley-Hamilton theorem 
[14 1 that for all i + j > n, one sets = J2h+k<n Ch,kA^’^’'^\ where Ch,kS. are real numbers. 
Therefore, for all i + j > n, C kerCA*^*’-^^ and consequently can be 

computed in a finite number of steps as, 

A^ = kerO= Pi (kerCA^*’^')). (20) 

*J>0, i+j<n 

Now, we consider the following subspace. 


^00 = 5 ^^ 


( 21 ) 


It follows that if xq = (• • • , x_i, Xo,Xi,- ■ ■) e then Xj G for all i e Z, and given the 
zero input assumption one gets y{k) = 0 for all fc G N (in equation @). By considering A^, 
where A is defined as in equation @ and A; G N, it can be shown that ^ A^g- Also, note 
that although is an Inf-D subspace, it can be computed in a finite number of steps (one 
only needs to compute A^). However, as explained in [Q, [ |M| the invariance property (this is 
addressed in the next subsection) of A^ is not lucid. Therefore, in the following a subspace of 
A^ is introduced such that it enjoys this geometric property. To define the subspace A^ one 
needs the following notation. 

Let us express A" to denote the sequence of multiplications of Ai and A 2 , where a is a 
multi-index parameter that specifies the sequence of the multiplication. For example, consider 
A“ = A 2 A 1 A 1 A 2 A 1 , where we have a = (2,1,1,2,1). The notation ||a|| denotes the number 
of all Ai and A 2 that are involved in the corresponding multiplication (for the above example, 
we have ||a|| = 5). Now, consider the following subspace (for more details on A^ refer to Q), 


A^= Pi kerC'A'’ 


( 22 ) 


a <n 


The following lemma shows that the subspace that is used in [30|, [361, [37 j as the unobserv¬ 
able (non-observable) subspace is indeed A^. 


Lemma 4. The subspace A^ can be computed in a finite number of steps according to the 
following algorithm, 


% = keiC % = ArVfe_i n A 2 n ker a 


(23) 
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Proof: First, note that Yi = ker C fl ^ ker CPiker C and P {AiA 2 )~^ ker C P 

( 742 ^ 1 )“^ kerCP kertPP kerC*. In other words, Yk = n||a||<fc(^°')~^ Note 

that for every pair of operators C : MA ^ W and F : — )■ M"^, one ean show that kertPF = 

F^^kertP. Therefore, it follows that Yi = A^s- This eompletes the proof of the lemma. ■ 
Now, we set P(^,oo = Note that although dim(iV 5 ^oo) = cxo, one ean eompute it in a 

finite number of step (by eomputing p(^). 

B. Ai^2-Invariant Subspaces 

As stated in the Subseetion |II-B[ the 2D system ([T]) ean be represented as an Inf-D system 
In order to formulate the eorresponding Inf-D invariant subspaees one needs the next two 
definitions. 


Definition 2. pJ]/ Consider the Inf-D system Q, where the operator A is bounded (according 
to Lemma^. The closed subspace %o G A’ = called A-invariant if AYoo ^ ■ 


Definition 3. [45] The subspace Y <zMA is said to be an Ai^ 2 -invariant subspace for the 2D 
system o if AiY -f A 2 Y C Y, where Ai and A 2 are defined as in equation 0- ■ 


Note that Y is Ai 2 -invariant if and only if it is invariant with respeet to Ai and A 2 (i.e. AfY C 
Y and A 2 Y C Y). The following theorem provides the eonneetion between the Definitions 
and[3l 


Theorem 3. Consider the 2D system 0 and the Inf-D system 0. Let Y^o = X] where 
Y C The subspace Y^o is A-invariant if and only if Y is Ai 2 -invariant. 


Proof: First, note that every x G Y^o can be expressed as x = 'fYk=-oo where x^ = 
■■■ , 0 , 0 ,a;^, 0 , 0 --- G Y^o and Xk^Y. Therefore, one only needs to show the result for x^. 
(If part): Assume Y is Ai 2-invariant. Consider the Inf-D vector x|. It follows that A,x| = 
[• ■ ■ , 0 , (^20;^)"’" , (AiXkp^ 0 , 0 ■ ■ ■ j""". Since Y is Ai^2-invariant, it follows that Ax^ G Y^o- 
(Only if part): Let fCYoo and Xq g Consecjuently, xq G . Since ^Axq — [* * * 5 ^7 

(^2X0))"*", (AiXo)"*", 0 , 0 , • ■ ■ j"*" G Yoo, it follows that AiXo G Y and A2X0 G Y, and consequently 
Y is Ai 2-invariant. This completes the proof of the theorem. ■ 

Consider the subspaces Y C MA and Y’ C If Y is the largest Ai_2-invariant subspace 
that is contained in 'Y, we denote Y =< 'Y\Ai^2 >■ We have shown in 0 that C A^, 
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and it is the largest Ai 2 -invariant subspaee that is eontained in ker C. Therefore, one ean write 
=< kerC'|74i^2 >■ Sinee ^ is 74i^2-invariant, by Theorem]^ ^,oo is ^-invariant. Therefore, 
if x(0) = (• • • , X-i, xo, xi, • • •) e ^s,oo (that is, x* G ^ for all i G Z) and zero input, 
x(/c) G ^,oo for alH G Z and y{k) = 0 for all A; G N (in equation (|^). We designate ^s,oo as 
the invariant unobservable subspaee. 


Remark 6 . As stated in Subsection II-B\ there are two different types of initial condition 
formulations. In this paper, we use the first formulation that is compatible with the Inf-D system 
©. Recall that the second formulation is expressed as x(i,0) = hiif) and x(0,j) = h2{i), 
where i,j G N. Now, let x(z,0) G ^ and x(0,j) G The Ai 2 -invariance property of NTg 
verifies that y{i,j) = 0. In other words, is also the invariant unobservable subspace 

of system o with the second formulation of the initial conditions. Therefore, without loss of 
any generality, one can apply our proposed approach to both initial condition formulations as 

Moreover, 1^ ihe largest A-invariant in the form A^S,00 = E ^ 


II-A 


provided in Section 
that is contained in ker C, where C is defined in 0 - 


C. Conditioned Invariant Subspaces 

Another important subspaee in the geometrie FDI toolbox is the eonditioned invariant (i.e., 
the (C, Ai 2 )-invariant) subspaee that is defined next. This definition is an extension of the one 


that has appeared and presented in |45| and [46|. 


Definition 4. The subspace #^3 = E ^ (where W C M”j is said to be the conditioned invariant 
subspace for the 2D system 0 if there exist two output injection maps Di, D 2 ^ M” such 
that {Ai-\- DiC)W C W and {A 2 -\- D 2 C)W C W. In other words, W is [A-\- DC]i^ 2 -lnvariant 
(i.e., invariant with respect to Ai+DiC and A 2 +D 2 C). We designate W as the finite conditioned 
invariant subspace (since dim(l^) < 00 ) of the 2D system Q. ■ 


Similar to ID systems, one ean now state the following result. 


Lemma 5. The following statements are equivalents. 

(i) The subspace is conditioned invariant. 

(ii) Ai(^nkerC') + A 2 (linkerC) C 
(iii) A(W^ n ker C) C Woo. 
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where ^ W. 


Proof: {i) {ii) and (i) {iu): By definition, there exists two maps Di and D 2 sueh 

that W is [A + DC] 1 , 2 -invariant. By utilizing Theorem is .Ad-invariant, where 


Ad = 


0 Ai + DiC A 2 + D 2 C 0 

0 0 Ai + D\C A 2 -\- D 2 C 


(24) 


By following along the same lines as in Lemma one ean show that Ad is bounded. Con¬ 
sequently, the result of ID system is also valid for the Inf-D system ([^. Henee, we have 
Aiifoo n kerC) C ifca (that shows (i) (Hi)). By eonsidering the strueture of A and C it 

follows that Ai{W fl ker C) + A 2 fW fl ker C) C if. 

{ill) ^ (i): Sinee A is bounded, the domain of A is equal to fb = therefore, the 

result of ID system is also valid for the Inf-D system @. Therefore, there exists a bounded 
operator V sueh that is A. -f DC-invariant. By eonsidering the strueture of and ker C, 
it is easy to show that one solution for D is given by 


V = 


0 Di Da 0 
0 0 Di Da 


(25) 


IS 


Henee, by using Theorem]^ the subspaee if is \Ai- DC] 1 , 2 -invariant and eonsequently 
a eonditioned invariant subspaee. This eompletes the proof of the lemma. ■ 

In the geometrie FDI approaeh, one is interested in eonditioned invariant subspaees that are 


eontaining a given subspaee [22|. Let us define all the eonditioned invariant subspaees eontaining 
a subspace (=5f C M”) as = \if^\ 3 Di and Da ; {Ai + DiC)if C 

if and if D .if}. It can be shown that for a given subspace (or the set n(-Sf) is closed 
under intersection, and hence the set 0.{Af) has a minimal member as iff^{^). The minimal 
conditioned invariant subspace containing a given subspace =^00 = (that is, if^fif)) is 

obtained by invoking the following non-decreasing algorithm that is provided below, 

if^ = Af, 

if'^ = S£ f (Ai(ir^-^ n ker C) + A2{if^-^ n ker C)), 


(26) 


September 20, 2016 


DRAFT 







22 


and W*, where W* = ko < n. 

Note that the above algorithm eonverges in a finite number of steps. Also, let be a finite 
eonditioned invariant subspace. The set of all maps D = [Di, Dq] such that W is [A + DC] 1 , 2 - 
invariant is designated by D{W). 


D. Unobservability Subspace 

The unobservability subspace p0| , [24| is the cornerstone of geometric FDI approach in ID 
systems. The following definition generalizes and extends this concept to the FMII 2D models. 
The Fin-D notion of this definition, for the first time in the literature, was introduced and utilized 
in Q for the FDI problem of 2D Roesser systems. 


Definition 5. The subspace Soo is said to be an unobservability subspace for the 2D system 0 
if there are three maps Di, D 2 , H such that S =< ker HC\[A + DC]i ,2 > and Soo = J^S. VFe 
designate S as the finite unobservability subspace of the 2D system Q. ■ 


Note that Soo is also conditioned invariant subspace and a generic unobservable subspace 
of the system (HC, Ai + DiC, A 2 + D 2 C). For accomplishing the gaol of the FDI task, one 


first computes an unobservability subspace and then obtains the matrix H [20|. Therefore, it 
is necessary to compute the unobservability subspace without having any knowledge of H. By 
following along the same lines as those used in 0 (Section IV.D), one can show that the limit 
of the following algorithm is the smallest unobservability subspace S* (and consequently S^) 
that contains a given subspace that is. 


= M*" and ^’^ = W* + n AY^^~^ n ker C), (27) 

where W* is the minimal finite conditioned invariant subspace containing .if, 5^ = ^ 5* and 

Finally, we set S^ = where S* = Finally, it should be noted that since W* + 

ker (7 = 5*+ ker C (this follows by applying ^”+1 = = S* in the algorithm ([TT])), one 

obtains D{W*) ^ D{S*). 

In this section, we first defined the invariance property of the 2D system 0 that are Inf-D 
subspace of the Inf-D system 0. Next, an invariance unobservable subspace (that is generically 
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equivalent to an unobservable subspace) ^,oo = X] ^ was introduced. Moreover, the con¬ 
ditioned and unobservability subspaces that are crucial in determining the solution to the FDI 
problem have been introduced. By utilizing the above results necessary and sufficient conditions 
for solvability of the FDI problem are subsequently derived and provided. 


IV. Necessary and Sufeicient Conditions eor Solvability of the FDI problem 


In this section, we first present sufficient conditions for detectability and isolability of faults. 
Next, by employing three different filters (namely, ordinary, delayed deadbeat, and LMI-based 
filters) sufficient conditions for solvability of the FDI problem are presented. 

Consider the faulty FMII model ([T]) (i.e., the system is subjected to two faults /i and / 2 ) and 
the detection filter @ designed to detect and isolate the fault /i. By augmenting the detection 
filter dynamics ([^ with the faulty 2D model ([T]), one obtains, 

Xe{i -l- 1, j -l- 1) = AiXe{i,j -|- 1) + A^Xeii + 1, j) “t" -|- 1) + -|- 1, j) + + 1) 

+ ^2,eh{'i + 1, j) + ^l,e/l(Li + 1) + -^l,e/l(* + Ij j): 


= CeXe{i,j), 


where Xe 
and. 




(o refers to the dimension of cu), Cg 



1 

o 

1_ 


1 

o 

c- 

1_ 


Bi 

Al = 

EiC Fi 

— 

? ^2 — 

E2C E2 

,Bt = 

Ki 


Bl = 


B2 

K2 


>Lie = 


(L: 


3\T 


, e {1,2}. 


HiC Ml 


(29) 


In this section, by assuming that the unobservable subspace of the above augmented system is 
Al 2 -invariant, it is shown that the sufficient condition is also necessary for solvability of the FDI 
problem. Moreover, an analytical comparison between our proposed approach and the method 


developed in p6| is also provided to highlight the strength and capabilities of our proposed 
methodology when compared to the literature. 


A. Main Results 

The following lemma provides an important property for the invariant unobservable subspace 


<00 (and that is associated with the system ([2^. 
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Lemma 6. Consider the 2D system ( [28| ) and its invariant unobservable subspace Then 

Q~^'y^s^oa unobservability subspace of the 2D system ©. with Q = diag(- ■ ■ ,Q,Q,- ■ ■) 
and Q representing the embedding operator into (i.e., Q : M" —>■ and Qx = 

Proof: First, recall that Note that, = ,^ = {x| [g] G CCf}, and 

assume that [o ] ^ According to the fact that is A® 2 -invariant 0, we have A® [g] G 
and if X G ker C then Aix G y. It follows that Ai(A^ Piker C) C By following along the same 
lines as those above one can show that A 2 (A^nker C) C Therefore, A^oo = Q~^^s% = Zl 
is a conditioned invariant subspace. Moreover, given that jCf is contained in kerCe, we have 
A^ C kerifiC. Therefore, the subspace A^ is a finite conditioned invariant subspace contained 
in keri^iC'. Since is the largest A® 2 -invariant subspace in kerCe, it follows that A^ is the 
largest [A + DC]i ^2 invariant subspace in keri^iC (i.e., A^ is a finite unobservability subspace 
of the 2D system ([T])). In other words, A^ is a finite unobservability subspace of the system ([T]). 
This completes the proof of the lemma. ■ 

The following theorem provides a single necessary and sufficient condition for detectability 
and isolability of faults (i.e., the existence of a subsystem such that it is decoupled from all 


faults but one - refer to Subsection II-C for more details). 


Theorem 4. Consider the 2D system o that is subject to two faults fi and / 2 . The fault fi 
is detectable and isolable (in the sense of Remark ^ if and only if the following condition is 
satisfied, 


(A^/ + ^ 5 * (30) 

where SI is the smallest finite uobservability subspace (refer to Definition of the 2D system 
0 containing + Af 2 ^ (this represents the limit of the algorithm ( [TT] ), where one sets Af = 
^2 "f -^2 the algorithm @). ■ 

Proof: (If part): Let Hi, D\ and D\ be the corresponding operators to (i.e. =< 
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keiHiC\A + DC >). Also, consider the following deteetion filter, governed aeeording to 

0Ji{i + 1, j + 1) =FiOJi{iC + 1) + F20Ji{i + 1, j) 

+ PiBiu{i,j + 1) + PiB2{i + 1, j), 

ri(i, j) + 1, j) - Hiy{iC), (31) 

where Pi : M"' —)■ iSj" is the eanonieal projeetion of M” on M"'/5i, Fk = + Do,kMi, 

k = 1,2 where AfPi = Pi(Ai + PiC) and A^Pi = Pi{A2 + D2C), Do^k are filter gains and 
Ml is the unique solution to Mi Pi = HiC. Now, by defining e{i,j) = Pix{i,j) — u:{i,j), one 
obtains 


e(* + 1, j + 1) = Fie{i,j + 1) + P 2 e(i + 1, j) 

+ Pifi{id) + Pifi{id), 

ri(i,j) = Mie(i,j). (32) 

Note that P 1 L 2 = P 1 L 2 = 0 and by the eondition ( [30| ), at least for /c = 1 and/or k = 2, we have 
PiL^ 7 ^ 0. Therefore, the residual signal ri is deeoupled from /2 and is sensitive to /i. In other 
words, fi is deteetable and isolable. 

(Only if part): We show the neeessary part by eontradietion. Let /i be deteetable and isolable, 
and C 5i. Henee, the residual signal ri is deeoupled from /2 and sensitive to fi. 

However, from Lemma one obtains + ^i) C Consequently, the fault fi is not 

deteetable by using the residual ri in ( [32l ), whieh is in eontradietion with the solvability of the 
FDI problem. ■ 

As stated in Remark Q the FDI problem has two main steps. Theorem [^provides the neeessary 
and suffieient eondition for the first step (that is, deteetability and isolability of /i). Therefore, 


eondition ( |3()| ) is also neeessary for solvability of the FDI problem. For the seeond step, we need 
to determine 'Do,k, and eonsequently F^, k = 1,2 in the filter pT]) sueh that the residual dynamies 


( |32| ) is asymptotieally stable. In what follows, we provide the eonditions that are based on three 
different observers (namely, the ordinary, the delayed deadbeat, and the LMI-based filters). We 
first need to present the following theorem. However, one first needs the following notation. 
Consider the PBH matrix ( [T^ . The matrix N{zi,Z 2 ) = [A^i(^i,^ 2 ), A^ 2 (^i,^ 2 )] is the minimal 
left annihilator of the PBH matrix if N{zi, Z 2 )PBH{zi, Z 2 ) = 0, N{z\, Z 2 ) is full row-rank, and 
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any other left annihilator Ni{zi, Z2) can be expressed as Ni = N^N, where is a polynomial 
matrix of zi and Z2- 


Theorem 5. / [77| / Consider the 2D system ©• The FDI problem is solvable (by using the 
approach that is proposed in ) if and only if, 

ZiL^ + Z2L‘^ 


iajik{N(zi, Z 2 )] 


Jqxp 


) = p , \/zi,Z2 e C - {0} 


(33) 


where p is the number of faults, N(zi,Z 2 ) = 


the PBH matrix ([TO]), and U = 


... 


iVi, N2 


denotes the minimal left annihilator of 


for k = 1,2. 


The following corollary provides the sufficient conditions for solvability of the FDI problem 
by utilizing our proposed method as well as a (delayed) deadbeat filter. 


Corollary 2. Consider the 2D system O that is subject to two faults fi and / 2 . Provided 
condition ( |30| ) is satisfied and condition of Theorem [^ is satisfied for the quotient subsystem 
(HiC,Af,A^), where Hi,A^ and A 2 are defined in then the FDI problem is solvable. 

Proof: By invoking Theorem]^ the residual signal ri in (jSTj) is decoupled from all faults but 
/i. Now, given that the condition ( [ 3 ^ is satisfied for the subsystem (HiC,A\,A^, by determining 
the observer gains Do,i and Vog, one can design a delayed deadbeat observer such that rfii — 
ni,j — 712) = (refer to [ [T 7 | Theorem 1 ]). Hence, the FDI problem is solvable. This 

complete the proof of the corollary. ■ 

The following corollary provides a sufficient condition for solvability of the FDI problem by 
utilizing an ordinary (without delay) deadbeat detection observer. 


Corollary 3. Consider the 2D system 0 that is subject to two faults fi and / 2 . Let the condition 
(jig) be satisfied and the PBH matrix of the quotient subsystem (HiC,Aj(,A^) be right zero prime, 
where Hi , Af and A 2 are defined in igg. Consequently, the FDI problem is solvable by using 
a deadbeat (without a delay) observer. 


Proof: By invoking Theorem]^ the residual signal in ( | 3 ^ is only affected by fi. Therefore, 
the detection and isolation of fi is reduced to that of designing Do 1 and Vog gains. Since the 
PBH matrix of the subsystem (HiC,AfiAfi) is right zero prime, by invoking Theorem one 
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can determine the gains sueh that the error dynamies ( |3^ is asymptotieally stable and the FDI 
problem is solvable. This eomplete the proof of the eorollary. ■ 

As shown in Theorem under eertain eonditions one ean obtain a residual signal that is 
deeoupled from all faults but one. Therefore, design of a residual generator to deteet and isolate 
the fault fi in the 2D system ([T]) is redueed to that of deteeting this fault in the system ( [3T] ) by 
using an observer. The Corollaries and [^provide suffieient eonditions for solvability of the FDI 
problem by using delayed and ordinary deadbeat filters, respeetively. However, as pointed out 
in [ [3^ , design of an observer for FMII models is based on polynomial matriees. This method is 
unfortunately not always numerieally or analytieally straightforward to develop and therefore, in 
this work we also develop another set of suffieient eonditions for solvability of the FDI problem 
by using a 2D Luenberger observer. 

The next eorollary provides suffieient conditions for solvability of the FDI problem by using 
a 2D Luenberger observer. 


Corollary 4. Consider the 2D model o. where the condition ( [30| ) is satisfied. The FDI problem 
is solvable if there exist two symmetric positive definite matrices Ri and R 2 such that, 

w2( ° - + ( 34 ) 

_o fi,] [(/i5)T 

where A\ and are defined in equation and the columns of Wm are the basis of ker M. 


Proof: 

By invoking Theorem]^ the fault /i is deteetable and isolable, and eonsequently the residual 
signal Ti in ( [3T] ) exists (and is deeoupled from all faults but /i). By eonsidering the LMI eondition 
( |34l ) and invoking results from Corollary the error dynamies ( |3^ is asymptotically stable. If 
/i(-, ■) = 0, the residual signal rfi-, •) eonverges to zero as z + j — )■ cxd. Otherwise, the residual 


has a value that is different from zero. Therefore, the eondition ( [7a| ) is satisfied and the FDI 
problem is solvable. By following along the same proeedure as those above one ean also design 
another state estimation observer to deteet and isolate the fault /2(-, •)• Therefore, this eompletes 
the proof of the eorollary. ■ 


Remark 7. It is worth noting that one can directly work with {as defined in ( [T ^ ) and derive 
necessary and sufficient conditions by following along the same steps as those that have been 
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proposed in / |?7| /. However, there are two main drawbacks associated with this approach that 
are as follows: 

1) The invariant subspaces are not necessarily computed in a finite number of steps. 

2) By factoring out .yVg, the resulting subsystem does not necessarily have a Fin-D represen¬ 


tation. For more clarification on 2D realization, refer to the work in [ 48 ]. 


To summarize, Theorem provides a single necessary sufficient condition for detectability 
and isolability of the faults (refer to Subsection |IVj) that is also necessary for solvability of the 
FDI problem. Three sets of sufficient conditions for solvability of the FDI problem by utilizing 
(i) a (delayed) deadbeat observer, (ii) an ordinary deadbeat observer, and (ii) a 2D Luenberger 
observer are derived in Corollaries and respectively. 

Table summarizes the main results that are developed and presented in this subsection. 


TABLE I: Pseudo-algorithm to detect and isolate the fault fi in the 2D system ([T]). 


1) Compute the minimal conditioned invariant subspace 'W* * containing all .SCj + 

subspaces such that j 4^ i (by invoking the algorithm where .5^ = + 

^!)). 

2) Compute the unobservability subspace S* containing (by using the 

algorithm <27}) . 

3) Compute the operators Di, D 2 such that W* is the minimal conditioned invariant 
subspace of the 2D system Q- 

4) Find the operator H such that keriifC' = 5* -t-kerC. 

5) If S* r\ 7 ^ 0 or S* C\ 7 ^ 0, then the fault /i is detectable and isolable (refer 
to Remark and Theorem and 

• If the conditions of the Corollary are satisfied, the FDI problem is 
solvable by using a (delayed) deadbeat filter. 

• If the conditions of the Corollary are satisfied, the FDI objective can 
be accomplished by using a deadbeat (without a delay) observer. 

• If the conditions of Corollary are satisfied, there exist an LMI-based 
observer for detection and isolation of the fault fi. 

6 ) The output norm of any of the above detection filters is the residual that 
satisfies the condition ( 0 . 
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B. Comparisons with the Other Available Approaches in the Literature 

In this subsection, our proposed approach is now compared and evaluated with respect to the 


existing geometric methods in the literature [361, [37|. We first show that if the FDI problem is 
solvable by using the approach in the above literature, our approach can also detect and isolate 
the faults. Furthermore, we provide a numerical example where it is shown that our approach 
is capable of detecting and isolating a fault, however, the necessary conditions provided in the 
algebraic methods in [ [T5| , [ fT7| as well as [ |^ are not satisfied. 

The equivalent 2D version of the necessary condition (as derived in [ [36| , Theorems 2 and 3) 
to detect and isolate the fault /i can be summarized as follows: The fault fi in the 2D system ([T]) 
is detectable and isolable according to if CW^* = 0, where and denote the 

minimal finite conditioned invariant subspace that contain + =2^^ and respectively 


(refer to [36|, Theorem 2). It should be pointed out that the observability assumption of (C, 


Ai, A 2 ) is a fundamental requirement and condition in [361 (although it was stated in [36j that 
this assumption was made for simplicity of the presentation). The main reason for the above 


limitation lies on and is due to the fact that the proposed approach in [36| is based on the results 


of [211. However, as stated in [21 [ the observability assumption is quite a crucial and critical 
condition. For further illustration and clarification for the above serious concerns consider the 
following 2D system. 


x{i + 1, j + 1) = 


y{h3) = 


0 1 
0 0 

0 1 


x{i,j + 1 ) + 


0 0 
0 1 


x{i + 1, j) + 


/i + 


/ 2 , 


(35) 


We have = kerC, W* = ^ 2 , and consequently C#)* fl CW 2 = 0. Therefore, 

the sufficient condition for solvability of the FDI problem under the zero initial condition in 
[ |3^ (Theorem 2) is satisfied. However, it is easy to verify that C JC, and consequently 
fl is not even detectable (in other words, fi has no effect on the output signal). Our proposed 
methodology does not suffer from the above limitation and restriction. 

We are now in a position to state the following theorem. 

Theorem 6. Consider the 2D system ([T]) and assume that the FDI problem is solvable by using 
the approach that is proposed in /[36]/. Then the approach that we have proposed in this work 
can also detect and isolate the faults in the system Q. 
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Proof: According to Theorem 3 in p^ , CWf fl CWf = 0, and Wf and Wf are inter¬ 
nally/externally stabilizable. Therefore, there exist two maps Di and D 2 sueh that and Wf are 
both [A+DC] 1 , 2 -invariant, and the pair {Ai+DiC, A 2 +D 2 C) is stable (that is, the eorresponding 
2D system is asymptotieally stable). Sinee CWf fl CWf = 0, there exists a map H sueh that 
HiCWf = CWf and HiCWf = 0. It follows that fl = 0, where is the invariant 
unobservable subspaee of {HiC, Ai + DiC, A 2 + D 2 C). Note that is an unobservability 
subspaee of (C, Ai, A 2 ) eontaining since is the smallest unobservability 

subspaee eontaining it follows that fl + Af^) = 0. Also, sinee (Ai + DiC, 

A 2 + D 2 C) is stable, it ean be shown that (A^, Af) is also stable, where P[Ai + DiC] = A\P 
and P is the eanonieal projeetion of Therefore, one ean eonstruet an observer for the residual 
system ( [3T] ) to deteet and isolate the fault fi (by ehoosing Doi = D 02 = 0). By following along 
the same proeedure, one ean deteet and isolate the fault / 2 . This completes the proof of the 
theorem. ■ 


Remark 8. Theorem shows that our proposed approach can detect and isolate faults that 
are detectable and isolable by using the geometric method in ( |3d| /. However, as shown below 
an example is provided where this approach fails whereas our proposed approach can still 
detect and isolate the faults. Also, we show that in this example the methods in / f75| /, / [77| / 
are also not applicable. 


1) Illustrative Example (Limitations of the Methods in /[75|/, / [TZ] / and /[3d|/).' 

Consider the 2D system ([T]) that is subjeeted to two faults fi and / 2 . As stated earlier, in this 
work our main eoneem is on aetuator faults and it is assumed that the output signals are not 
affeeted by faults (refer to Remark [^. 

Consider the 2D system ([T]) that is subjeeted to two faults fi and /2 where, 

02x2 0 

I I 


Ai — 


00 n cr 
0 0.5 


0 0 


2 x2 


, A 2 = 0.5 


, L\ — L 2 — 0 , Bi — B 2 — 0 


L} = [0,0,0,1]T, LH [0,0,-1,1]^, C' = 


10 0 0 
0 0 0 1 


(36) 
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The FDI problem is solvable by using the approaeh in p3| | if and only if, 

- ziAi — Z 2 A 2 zi[L\, Ll,] 


rank( 


) = n , 'izi,Z 2 e C 


(37) 


C 0 

The full rankness of the PBH matrix ( [T^ is the neeessary eondition for the above property p3| 


(page 231-item (i)). However, for the 2D system (36), one gets rank(Pi?i7(l, 0)) = 3 < 4. 


Therefore, the FDI problem is not solvable by using the method in [15|. 

To eheek the eonditions of Theorem let us first obtain N{zi,Z 2 ). Consider [a, 6, c, d, e,/] 
as a row of N{zi,Z 2 ). Since N{zi, Z 2 )PBH{zi, Z 2 ) = 0, it follows that. 


0.52:16 + (1 — 0 . 5 z 2 )d + / = 0 
{l-0.5zi)b-0.5z2d = 0 


(38) 


a — 0.5 z 2C + e = 0 

— 0.52:10 + (1 — 0 . 5 ^ 2 )c = 0 
From the first set of equations, one choice for a, c and e is given by a = 2 — 2 : 2 , c = 2:1 
and e = 0 . 52 : 2^1 + Z 2 — 2, respectively. Also, based on the second set one can write b = Z 2 , 


d = 2 — 2:1 and / = 2 — 2:2 —; 2 i. It can be shown that A^( 2 :i, 2 : 2 ) = 


a 0 c 0 e 0 

0 b 0 d 0 f 


Therefore, 


N{zi,Z2) 


Zi[Lm] 


0 


= 


U 2:1 

2 — 2:1 2-2:1 


. It is easy to check that the rank(A^( 2 : 1 , 2) 


Zi[L\,Ll] 

0 


2 for all 2:2 G C and 2:1 = 2. Hence, the necessary condition that is proposed in [17| (as given 


by equation (|3^) is not satisfied . 

Furthermore, the necessary condition to detect and isolate the fault fi by using the approach in 
|36| is (7#)* n CW 2 = 0. Since ^ ker(7, by invoking the algorithm ( |26l ), one obtains 

Wi = A£x and = ^£ 2 - It follows that, D = span{[0, I]"""} . Therefore, the 


necessary condition in [36| is also not satisfied. In other words, the fault /i cannot be detected 
and isolated by using the detection filter ([^, if one restricts the filter to the case with M = C 


(or H = I), according to the required results in [36|. 

Finally, it is now shown and demonstrated that one can detect and isolate both faults /i and 
/2 by using our proposed methodology. Towards this end, by invoking the algorithm ( [27] ), one 
can write iSj" = (that is, the finite unobservability subspace containing Af^) that satisfies 


the condition (30). By considering Di = 


-I T 


0 

0.5 


0 0 0 
-0.5 0 0 


and Do = 


n T 


0 0 0 
0 0 0.5 


0 

-0.5 


is [A -f PC] 1 , 2 -invariant. Also, since kerC -f 5]" = kerPiC and MP 2 = HC, one gets 


Hi = [1, 0] and Mi = [1,0,0]. Hence, the residual ( |32l ) that is only affected by the fault fi 
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'o 

0 

72' 

2 


’ 0 

0 

0 ' 

is given by A\ = 

0 

0.5 

0 


0 

0 

0 


0 

0 

0 


V2 

. 2 

V2 

2 

0.5 


and Bi = B 2 = 0. Since the pair (A^, 

A\) is stable, by eonsidering Doi = -D 02 = 0, the deteetion filter for the fault fi (as given by 
equation @) is now obtained aeeording to, 

+ 1, J + 1) = + 1) + + 1, j), 

(39) 

By following along the same procedure, one ean also design a deteetion filter to deteet and 
isolate the fault / 2 . Therefore, our proposed approaeh ean aeeomplish the FDI objeetives while 
the approaehes that are proposed in [fTSl, [[T7i| and [[3^ cannot achieve this goal. 

Remark 9. All the conditions for solvability of the FDI problem in the literature (for both ID 
and 2D systems) and also our proposed conditions are generic, although this fact is not explicitly 
mentioned. For clarification, consider the faulty model where x G Ai = A 2 = 0.4 * I, 
B\ = L\ = [1,1]"*", B 2 = L 2 = [0,1] and C = [1, —1]"'". Let the initial condition x(0) = 0 
and fi{i,j) = 1 for all i + j > 0. It follows that y = 0, and consequently fi is not detectable. 
However, below we show that sufficient conditions in the literature [15], as well as 

our proposed conditions are still all satisfied. 

1) (Conditions in ffT^ and 1151). It follows that N{zi, Z 2 )PHB{zi, Z 2 ) = 0, where N{zi, Z 2 ) = 

"-1 0 0.4(z, + Z2)-1 , , , 7 . , 

IS also satisfied. 


, and consequently the condition in Theorem 

0 1 0.4(^i + Z 2 ) — 1 

2) ( Conditions in ). By following the algorithm ( [2^ we obtain Wf = tmd Wf = F^ 2 - 


It follows that CWf nCWf = 0, and consequently the condition in [36] is also satisfied. 


3) (Our proposed Conditions) By following the algorithm ([27|) we obtain = Ffi and 


S 2 = L follows that S] fl Jf 2 = 0, and consequently the condition ( [30| ) is also 
satisfied. 

To summarize, in this seetion we have developed and presented a solution to the FDI problem 
of 2D systems by invoking an Inf-D framework for the first time in the literature and by utilizing 
invariant subspaees and derived neeessary and suffieient eonditions for solvability of the problem. 
It was shown that if the suffieient conditions for solvability of the FDI problem that are provided 


in [361, [371 are satisfied, then our proposed approaeh ean also deteet and isolate the faults. 
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However, as shown above there are eertain systems that the methods in p3| , p7| , are not 
applieable and eapable of deteeting and isolating faults, whereas our proposed approaeh ean 
solve this problem sueeessfully. 


V. Simulation Results 


In this seetion, we apply our proposed FDI methodology to a two-line parallel heat exehanger 
system 0, [ |^ . Speeifieally, we verify the neeessary and suffieient eonditions that are derived 
in the previous seetion, and also design a set of filters to deteet and isolate faults under both 
full and partial state measurement seenarios (this is to be realized by an appropriate selection 
of the output matrix C). 

The heat exchanger is usually subject to two different types of faults, namely the fouling 
and the leakage [ |49| . The mathematical model of a typical heat exchanger is governed by the 
following hyperbolic PDEs, 


dTf _ ^ dTf 


dt 


- -«/ 


dz 




BT BT 


(40) 


where Tj and Tg denote the temperature of the cold (fuel) and the hot (exhaust gas) sections, 
respectively. The coefficients «/ and ag are proportionally dependent on the speed of the fluid 
and the gas, respectively, and the coefficient (3 is related to the heat transfer coefficient of the 


wall [ |49| . Moreover, /i and /2 denote the leakage and the fouling effects, respectively. Finally, 
it is assumed that the boundary conditions (the inlet temperature) Tf(t, 0) and Tg{t, 0) are given 
and only the outer section (i.e. Tg) is subject to the leakage. 


A. The Approximation of Hyperbolic PDE Systems by 2D Models 

Let us first illustrate and demonstrate how one can approximate a general hyperbolic PDE 

system by using 2D models and representations. Consider the following hyperbolic PDE system 

Br^ ~ Br* ~ ~ ^ ^ ^ 

T A 2 X + Bu+ Lkfk, (41) 

fc=i 

where z denotes the spatial coordinate, x{z,t) G u{z,t) G and fkiz,t) ^ denote the 
state, input and fault signals, respectively. Also, the operators Ai, A 2 , B and Lfs are real matrices 
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with appropriate dimensions. Note that every hyperbolie PDE system with eonstant eoeffieients 


ean be represented in the form ( |4T] ) with a diagonalizable Ai [50| (Chapter 1, Detention 1.1.1). 
By applying the finite differenee method to the system ( [4T| ), one obtains 

5:(iAz, (j + l)At) — i(zAz, jAi) ~ x{iAz,jAt) — x{{i—l)Az,jAt) 


At 


—Ai- 


Az 


■ A2x(iAz, j At) 


Bu{iAz,jAt) + y^Lkfk{iAz,jAt). 


(42) 


fc=i 


By setting x{i,j) = 


x{{i — l)Az,jAt) 
x{iAz,jAt) 


, we ean now write 


X[ 


k=l 


'{i + Ij + 1) = Ax{i,j + 1) + A2x{i + 1, j) + Biu{i + l,j) + ^Lkfk{i + 1, j), (43) 

where u{i + 1, j) = u{i,j) for all i and j, Ai = 


0 / 
0 0 


and 



1 

0 

0 

_1 


0 


0 

A2 — 

(/+ If + Ati2)_ 

, B,= 

B 

, Lk — 

-1 


, fihj) = + l,j)- 

Therefore, the PDE system ([4T]) ean be approximated by the EMU 2D model ([4^. 


Remark 10. As shown in [10], the parabolic PDE systems can be approximated by FMII 3D 


models. As shown in Appendix^ our proposed EDI methodology for 2D systems can be extended 
to 3D systems, therefore the extension of the results of our proposed EDI methodology can also 
be applied to parabolic PDE systems. ■ 


Eor the purpose of eondueting simulations, the parameters in equation (|40|) are taken as a/ = 


ag = [3 = 1. Also, by eonsidering Az = At = 0.1 one ean diseretize the system (40) as. 


fi + l,j + 1) = 


O 2 


x2 J-2X2 


x{i,j + 1) + 


02> 


02> 


r - 0.1 0.1 

12 x 2 0.1 - 0.1 


fi + 1, j) + Lifif + l,j) 


(44) 


02x2 02x2 

+ B\u{i + 1, j) + L\f 2 {i + 1, j), 
yfA) =Cx{i,j), 

where Lj = [0, 0, 0, l]"^, = [0, 0,-1,1]'^ and x(i, j) = [Tg{{i-l)Az,jAt),Tf{{i-l)AzjAt), 

Tg{iAz,jAt),Tf{iAz,jAt)]^. Einally, by following along the same steps as in |j^, one obtains 

02x2 


Bo = 


1.1 -0.1 
-0.1 1.1 


, m( 0, j) = [T/(0, Af),Tg(0, Af)]^ and u{ij) = u{0,j) for all i. 
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We first assume that both temperatures Tf and Tg are available for measurement along the 
spatial eoordinates at diserete points (i.e., Tf{iAz,jAt) and Tg{iAz,jAt) are available from 
sensors). Next, we eonsider the ease where only the outer temperature (that is, Tg) is available 
for measurement. As we shall show subsequently, in the latter ease by using the ID approximate 


ODE model (for example, as in [511), the faults fi and /2 are not isolable, whereas by using 
our proposed 2D EDI methodology one ean deteet and isolate both faults. 


0 0 0 1 


B. FDI of a Heat Exchanger by Using Full State Measurements 

Eet us assume that both temperatures (namely, Tj and Tg ) are available for measurement. 

10 0 0 

Therefore, one ean seleet the output matrix as O' = 

1) Detectability and Isolability Conditions: 

By applying the algorithm ( [26l ), the minimal finite eonditioned invariant subspaee eontaining 
^ = spanjL^} is obtained as Wf = and by applying the algorithm ( [27] ), one obtains 
S* = Wf = ^ 2 - It follows that Tfl fl 5* = 0. Therefore, the suffieient eondition (as given by 
Theorem]^ is satisfied. By following along the same lines the neeessary and suffieient eondition 
for deteetability and isolability of the fault /2 ean also be shown to be satisfied. 

2) FDI 2D Luenberger Filter Design: 

As stated earlier, we are interested in designing 2D Euenberger deteetion filters by using the 


EMI eondition that is proposed in the Subseetion 11-E In this part of the paper, we design a filter 
for deteeting and isolating the fault /i (without loss of any generality, by following along the 
same lines as eondueted below one ean also design a filter to deteet and isolate the fault / 2 ). The 
2D deteetion filter must be deeoupled from the fault /2 (refer to the eonditions in equation Q). 
As stated above, the finite unobservability subspaees eontaining the subspaee ^2 obtained by 
using the algorithm ( |27| ) and is given by S* = Wf = 

The output injeetion matriees D\ and D 2 for Wf are to be derived sueh that Wf is [A + 
DC] 1 , 2 -invariant. Therefore, one ean write {Ai + DiC)W 2 = 0 and {A 2 + D 2 C)W 2 = 0, 
where the eolumns of W2 are the basis of Wf fl fWf fl ker C)-*“. One solution to and D2 is 


Di = 


0 0 0 0 
1-10 0 


and Do = 


0 0 0 0 

0 0 - 0.2 0.2 


. Also, let P = 


10 0 0 
0 10 0 
0 0 11 


is the eanonieal projeetor of the subspaee 5C), where P is used in equation ( |3T] ). By using 


(whieh 
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'kei HiC = + kerC and Mi Pi = HiC, we have H = [1, 0, 0] and the output matrix M 

beeomes M = [1,0,0]. Henee, the faetored out 2D system is now expressed as, 


uj{i + l,j + 1 ) = 


0 0 - 1.42 
0 0 0 
0 0 0 


+ 1 ) + 


0 0 0 
0 0 0 
- 0.7 - 0.7 0 


uj{i + l,j) + 


0 0 
0 0 
0.7 0 . 7 . 


u(* + i,j)+ 1 /i(i + i,j), 


ypi'iJ) = [i,o,o]w(k j), 

where uj{i,j) = Pix{i,j). It is straightforward to show that the positive definite matriees Pi = 
diag(0.4,1,2.133) and P 2 = diag(0.4,2.15,0.86) satisfy the inequality WjAcWc^ < 0. By 


using Remark one ean obtain Doi = 0 and D 02 = 
deteet and isolate the fault fi is given by, 


0 0 -0.7 


. Therefore, the filter to 


+ 1,J + 1) = 


0 0 - 1.42 
0 0 0 

0 0 0 


+ 1 ) + 

Do2P2y{i,j + 1 ), 


0 0 0 
0 0 0 
- 0.7 - 0.7 0 


+ l,j) + 


0 0 
0 0 
L0.7 0.7J 


u{i + 1, j) 


nihj) = [l,0,0]a;(z, j) - Pr/(z,j). 

Designing a filter to deteet and isolate the fault /2 follows along the same lines as those given 
above for the fault /i. These details are not ineluded for brevity. 

3) Threshold Computation: 

Due to presenee of input and output noise, disturbanees, and uneertainties in the model, the 
value of the residual rk{i,j) is not exaetly equal to zero under the fault free situation. Therefore, 
to reduee the number of false alarm flags one needs to apply threshold bands to the residual 
signals. In this subseetion, we present an approaeh for determining the thresholds that are needed 
for aehieving the FDI task. For ID systems, there are a number of approaehes for eomputing 
a threshold, e.g. based on the maximum or the root mean square (RMS) of the residual signal 
[ |5^ . In this work, we use the maximum residual norm. However, one ean also apply the RMS 
approaeh to 2D systems. 

Consider the FMII 2D model ([T]) subjeet to the fault free situation. The threshold th^ is then 
determined from, 

fhfc = Max \ri{i,j)\ , fori,j<iVi (45) 


where | ■ | denotes a norm funetion (in this work we use the norm-2), is the number of the 

denotes 


Monte Carlo simulations (refer to [53|) that are used to determine the threshold, r 
the signal of residual in the length and Ni is a suffleiently large number. 


k\ I 
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(a) The residual signal ri for detecting and 
isolating the fault fi. 



(b) The residual signal r2 for detecting and 
isolating the fault /2. 


Fig. 1: The residual signals for deteeting and isolating the fault /i for the first seenario (a single 
fault ease) using full state measurement. 


By utilizing predefined thresholds that are denoted by thk^k = 1,2, the FDI logic can be 
summarized as follows, 


if ri > thi ^ the fault /i has occurred. 


(46) 


if r 2 > th 2 ^ the fault /2 has occurred. 

Let us now consider two scenarios. In the first scenario, a single fault fi with the severity of 1 
occurs at (7 = 5 and j > 60). In the second scenario, multiple faults /i and /2 with the severity 
of 1 occur at (f = 5 and j > 50) and (i = 5 and j > 70), respectively. The residuals ri and r 2 for 
the first scenario are shown in Figure For the second scenario, the results are shown in Figure 
1^ The thresholds are determined by conducting Monte Carlo simulations [ (5^ corresponding to 
the healthy 2D system. As shown in Figures and our proposed methodology can detect and 
isolate the faults in both single- and multiple-fault scenarios (according to the FDI logic that is 
given by equation (|46l)). 


C. FDI of a Heat Exchanger Using Partial State Measurements 

In this section, we assume that only the outer temperature (Tg) is available for measurement. 
This corresponds to a more practical and physically feasible scenario in various applications 
(sensing the inner temperature requires a more sophisticated hardware). In this case, we set 
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(a) The residual signal ri for detecting and isolating (b) The residual signal r2 for detecting and isolating 
the fault /i. the fault /2. 

Fig. 2: The residual signals for deteeting and isolating the faults /i and /2 for the seeond seenario 
(a multiple fault ease) using full state measurement. 


C = 


(in other words, we measure Tg{i — 1, j) and Tg{i,j)). In this subseetion, 


0 10 0 
0 0 0 1 

we demonstrate and illustrate the eapabilities of our proposed FDI approaeh based on the 2D 


system modeling, whereas it is shown that by using a ID approximation of the PDF system ( |40| ) 
the FDI problem eannot be solved. 

1) ID Approximation of the Heat Exchanger System: 


The hyperbolie PDE system ( |4()| ) ean be approximated by applying diseretization through the 2 ; 
eoordinates as follows. Let i denote the length of the heat exehanger that is diseretized into N 
equal intervals (i.e., Az = -^). By using the approximation _ Tg{kAz,t)-Tg{(k-i)Az,t) 


dz 


Az 


one ean represent the PDE system (40) by the following ID approximate model. 


N 


X = 


Ax + Bu + + pg), 

y = Cx, 


k=l 


(47) 


Tg{Az, t),Tf{Az, f), • • ■ , Tg{NAz, t),Tf{NAz, f)] ^ e B = 

1,0,0,0,-- 

■ ,0 


0,1,0,0,-- 

■ , 0 _ 


where x{f) = 


and /f (f) = fi{kAz, t) (i = 1, 2). Also, Lf = [0 , • - , 0 , —1,1,0, • • ■ ]^, and the fault signatures 

(fc-i) 

L 2 are 2A-dimensional veetors sueh that only the element is 1 and the rest are zeros. 
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Moreover, 


yl = 



0 

0 ••• 











1 

0 

0 

0 


A 2 


0 ••• 

,C = 

0 

0 

1 

0 


0 

A 2 

A, ••• 

0 

0 

0 

0 



0 

0 









(48) 


in whieh Ai = 

It ean be shown that 


1+Az 

1 


1 

1 

0 

_1 

Az 

and A 2 = 

1 

1+Az 

Az 


0 -yf 

Az_ 


= span{L 5 ^,L 2 } , where 


. Now, eonsider the faults and /|. 
is the minimal eonditioned invariant 


subspaee eontaining (from the ID system perspeetive). Consequently, fi Afi ^ 0 
denotes the unobservability subspaee eontaining ^2 the ID system sense). Consequently, the 
faults fi and /| are not isolable. However, we show below that the faults ean be deteeted and 


isolated if one approximates the system ( |40| ) by using the 2D model representation. 

2 ) 2D Representation of the Heat Exchanger: 

Let us set and X 2 {iJ) = so that the system (|4^ 

is approximated by the system ( |4^ where all the operators are defined as before exeept for 


L\ = 


0 10 0 


and C = 


0 10 0 
0 0 0 1 


Note that sinee only the state Tg is assumed to be 


available for measurement, we sense Tg{i — 1, j) and Tg{i,j). By applying the algorithm pT]), 


where Af = one gets S* = spanjL^, 


0 0 11 


} . Sinee, = 0, 


10 0 0 

the fault /i is both deteetable and isolable (aeeording to Theorem Q. It should be noted that 
by applying the approaehes in p3| , p7| , p^ , one can also detect and isolate the fault /i. By 
following along the same lines as the ones given earlier one can show that the fault /2 is also 
detectable and isolable, where one can design the required detection filters. Figures]^ and [^depict 
the simulation results for the scenarios that presented above. As can be observed from Figures 
and 1^ (particularly. Figures and [4F]), it follows that fewer available information (recall that 
both Tf and Tg are measurable in Figure]^ whereas in Figure]^ only Tg is measurable) results 
in a situation where the spatial coordinate of a fault cannot be estimated accurately. 


VI. Conclusions 

The fault detection and isolation (FDI) problem for discrete-time two-dimensional (2D) sys¬ 
tems represented by the Fomasini-Marchesini model II (where other 2D models are subclasses 
of this representation) is investigated in this work. In order to derive the necessary and sufficient 
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(a) The residual signal ri for detecting and (b) The residual signal r2 for detecting and 
isolating the fault fi. isolating the fault /2. 


Fig. 3: The residual signals for deteeting and isolating the fault /i for the first seenario (a single 
fault ease) using partial state measurement. 




(a) The residual signal ri for detecting and isolating (b) The residual signal r2 for detecting and isolating 
the fault /i. the fault /2. 

Fig. 4: The residual signals for detecting and isolating the faults fi and /2 for the second scenario 
(a multiple faults) using partial state measurement. 


conditions for solvability of the FDI problem, the notion of the conditioned invariant and 
unobservability subspace of ID systems was generalized to 2D systems by using an infinite 
dimensional (Inf-D) framework and representation. Moreover, algorithms for computing and 
constructing these subspaces are introduced and provided that converge in a finite and known 
number of steps. By applying the linear matrix inequality (LMI) approach, sufficient conditions 
for existence of an asymptotically convergent 2D state estimation observer is derived. Necessary 
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and Sufficient conditions for solvability of the FDI problem are also provided. Furthermore, 
another set of sufficient conditions (that are based on deadbeat observers) are also provided. 
It was shown that although the sufficient conditions for applicability of the currently available 
geometric results in the literature are also sufficient for our proposed approach to accomplish 
the FDI goal, however, there are 2D systems where the geometric approaches in the literature 
are not applicable to detect and isolate the faults, whereas our approach can still achieve the FDI 
objective and task. Finally, simulation results are provided for the application of our proposed 
FDI methodology to an heat exchanger system to demonstrate and illustrate the capabilities and 
advantages of our proposed solution as compared to the alternative ID representation and ID 
FDI approaches. 


Appendix 


The results that are derived in this paper are dependent on invariant concepts that are provided 
in previous sections. Since these invariant subspaces are independent of the order of A1 and 
A2, there is a natural extension that is possible to FMII 3D models. In this appendix, we briefly 
show how our previous results can be extended to 3D systems. Furthermore, it should be pointed 
out that one can also extend these results to nD systems 
Consider the following 3D FMII model 

x{i + 1, j + 1, /c + 1) = Aix{i,j + 1, fc + 1) + A 2 x{i + 1, j, fc + 1) 

+ A 3 x{i + 1, j + 1, k) + Biu{i,j + 1, fc + 1) 

(49) 

+ i?2u(i + 1, J, k + 1) + B3u{i + 1, j + 1, k), 
y{i,j,k) = Cx{i,j,k). 

The corresponding Inf-D system is given by 


x(A; + 1) = Ax.{k), k eN 
y{k) = (3x(A;), 


(50) 


where x(A;) e Af = Xl 1^”. y{k) = (■■•, l/(-l + k, 1)'^, y{k, 0)'^, y{l + k, -1)"^, • • • )'^ e X 
and A is an Inf-D matrix with Ai, A 2 and A 3 as diagonal and upper diagonal blocks, respectively. 
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with the remaining elements set to zero, and C = diag(- ■ • , C, (7, • ■ ■In other words, we have. 


A = 

... 0 

Ai 

A2 

A3 

0 ■■■ 

,c = 

... 0 

c 

0 ■■■ 


... 0 

0 

Ai 

A2 

A3 ■ ■ ■ 

... 0 

0 

c ■■■ 


The generalization of the Definition ean be stated as follows. 

Definition 6. The subspace Y is called Ai^ 2 , 3 -invariant if AiY C Y for all f = 1,2,3. ■ 


Also, it ean be shown that the finite unobservable subspaee of the 3D system ( |49l ) is given 
by 

where 


^( 0 , 0 , 0 ) _ J = 0 if z or J or fc < 0 . 

In view of the above, one ean now extend the invariant unobservable subspaee to = 
ri||a||>oker CA", where, a is a multi-index parameter, and |a| and ||a|| are defined aeeording 
to the same manner that were defined for 2D systems (refer to the notation in the Introduetion 
seetion). For instanee, if a = (3, 1 , 3, 2, 1 , 1 , 1 , 3, 3, 3), then |a| = (4, 1 , 5) and | |a| | = 10 . Also, 
we have A^,oo = Y. Yds- 

By following along the same lines as in Theorem and Lemma we have the following 
eorollaries. 


Corollary 5. Consider the 3D system ( |4^ and the Inf-D system ( |50l ). Let Y^o = where 

Y C M". The subspace Y^a is A-invariant if and only if Y is A 12 , 3 -invariant. 


Corollary 6. Consider the 3D system ( |49l ) and the Inf-D system ( [50| ). The following statements 
are equivalents. 

(i) The subspace is conditioned invariant. 

(ii) Ai (irn ker C) + A 2 (ir n ker C) + A 3 ( ir n ker C) c W. 

(iii) AiWoo n ker C) C W^. 
where = Y ■ 
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By following along the same steps as above, the invariant subspace and the results provided 


in Sections III and can be extended to nD systems. Therefore, one can generalize the results 


of this paper to nD systems without major barriers and challenges. 
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